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Abstract. We study the properties of one-dimensional hypergeometric integral so- 
lutions of the (/-difference ( "quantum" ) analogue of the Knizhnik-Zamolodchikov- 
Bernard equations on tori. We show that they also obey a difference KZB heat equation 
in the modular parameter, give formulae for modular transformations, and prove a com- 
pleteness result, by showing that the associated Fourier transform is invertible. These 
results are based on SL(3,Z) transformation properties parallel to those of elliptic 
gamma functions. 



1. Introduction 



FTVl 



In this paper we continue the study of the q-analogue of the Knizhnik-Zamolodchikov- 
Bernard (qKZB) equations on elhptic curves and their solutions initiated in 

Fvl . 

In [FTVl], hypergeometric solutions of qKZB equations were introduced. In 



FTV2 



the monodromy of hypergeometric solutions was calculated, and a symmetry between 
equations and monodromies was discovered: the equations giving the monodromy are 
again qKZB equations but with modular parameter and step of the difference equations 
exchanged. In ||FV1|| we introduced the g-analogue of the KZB heat equation, which 
governs the change in the modular parameter of the elliptic curve, proved that it is 
compatible with the other equations and explained how to recover the KZB heat equation 
in the semi classical limit. 

In this paper we prove three results about our hypergeometric solutions in the case 
where the sum of highest weights is two: we show that the hypergeometric solutions 



also obeys the qKZB heat equation of |[F Vl|] , see Theorem We give a formula 
(Theorem |3.3|) for the transformation properties of the hypergeometric solutions under 
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the modular group SL(2, Z). We prove a completeness result, Corollary by showing 
that the associated "Fourier transform" is invertible. 

Then we show that these results are part of a bigger picture, in which the modular 
group combines with the transformation defined by the qKZB heat equation to give 
a set of quadratic identities for our generalized hyper geometric integrals. In fact, this 
picture can already be seen in a simpler situation, in which the elliptic gamma function 
0, ||FV2|| plays the role of the hypergeometric integral. The elliptic gamma function 



is a function T{z,T,p) of three complex variables obeying identities |[1''V2|| involving its 
values at points related by an action of SL(3,Z) k Z^. These identities mean that F is 
a "degree 1" generalization of a Jacobi modular function, see [ |1^'V2| ]. 

These identities are a scalar version of the identities obeyed by the hypergeometric 
solutions of the qKZB equation. Their meaning is that the hypergeometric integrals 
define a discrete projectively flat SL(3, Z)-connection (i.e., a lift of the action to the pro- 
jectivization) on a vector bundle over "regular" orbits of SL(3, Z) acting on the variety of 
pairs (point in C'^ — {0}, plane through containing it). This is the content of Theorem 
|6.8| . The results on the elliptic gamma function are used in the proof, since the "phase 
function" which appears in the integrand of hypergeometric solutions is a ratio of gamma 
functions. In fact, we see "experimentally" that there seems to be a principle stating 
that to every identity obeyed by the gamma function, there corresponds an identity for 
hypergeometric integrals. The proofs of the identities consist in applying the gamma 
function identity to the phase function in the integrand, and then use a version of Stokes' 
theorem to relate the integrals. The second step is relatively simple in the case of the 
one-dimensional integrals to which we restrict ourselves here, but becomes exceedingly 
involved in the higher dimensional case. Proving our identities in the higher dimensional 
case, i.e., if the sum of highest weights is larger than two, remains a challenging open 
problem. An alternative approach to this problem is based on representation theory: in 
||FV|| a representation-theoretic interpretation of a degenerate version of the qKZB equa- 
tions was established. It was shown that traces of intertwining operators for quantum 
groups satisfy a version of the qKZB equations and are eigenfunctions of analogues of 
Macdonald operators. This fact indicates that our Theorem |2.1| i s an elliptic analogue of 
the Macdonald- Mehta identities proved by Cherednik, Q, ||EK|| , and the present work 
is a glimpse into an elliptic version of Macdonald's theory. 

In fact Theorem |6]^ concerns the case of generic parameters, where an infinite- 
dimensional vector bundle is preserved by the projectively fiat connection. In a next 
paper we will restrict our attention to special SL(3, Z)-orbits. The projectively fiat 
connection can be then defined on a finite-dimensional vector bundle of theta-functions 
which are a g-deformed version of the space of conformal blocks in conformal field the- 
ory. In this setting the analogy with Macdonald's theory will appear more explicitly. 
Another degeneration of the SL(3, Z) symmetry of our hypergeometric integrals are the 
SL(3,Z) symmetries of the ordinary Fourier transform indicated in |[I''V4|| . 
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2. Hypergeometric solutions of the qKZB equations 

We use the definitions and notations of [^VT|. Tlie elliptic s/2 qKZB equations are 
a compatible system of difference equations for a function v{z, A, r) of £ G C", A G C 
and r, Imr > 0, taking values in the zero weight space Va_[0] of a tensor product of 
Er,rf{sl2) evaluation modules. This space comes with a basis of eigenvectors of commuting 
operators h^'^\ i = 1, . . . ,n and depends on parameters A = (Ai, . . . , A„) G C". The 
qKZB equations have the form 

(1) v{z + p5i,T) = Ki{z,T,p,r])v{z,T), i = l,...,n. 
They are supplemented by the qKZB heat equation 

(2) v{z,t) =T{z,T,p,ri)v{z,T + p). 

The step p is a complex parameter and (5j)j=i^.,, is the standard basis of C". Here v 
is viewed as a function of z and r with values in the space ^(Va[0]) of meromorphic 
VA[0]-valued functions of a complex variable A. The qKZB operators Ki{z,T,p,rj) are 
difference operators in A and can be expressed in terms of a product of (dynamical) 
i?-matrices. The last equation is the qKZB heat equation and involves the integral 
operator T{z,T,p,rj). The latter is expressed in terms of hypergeometric integrals: let 
u{z, X, fi,T,p,r]) G Va[0] (S> Va[0] be the universal hypergeometric function as in [[FVl 



We may view it as a function r, p, 77) taking its values in the Va[0] ® VA[0]-valued 
functions of the two "dynamical variables" A and /i. Then it is a projective solution 
(i.e., a solution up to a constant factor) of the qKZB equations in the first factor, and 
of the mirror qKZB equations in the second: 

u{z + 6jP,T,p,r]) = Kj{z,T,p,r]) (g) Dju{z,T,p,r]), 
(3) u{z + 6jT,T,p,ri) = ® Kj{z,p,T,r])u{z,T,p,ri), 

u{z + 6j,T,p,r]) = u{z,T,p,r]). 

The mirror qKZB operators Ry {z,p,T,ri) are obtained from the qKZB operators by 
exchanging r and p and "reversing the order of factors", see |^Vl]. The operators Dj, 
are operators of multiplication by certain functions of the h^^^ and the dynamical 
variable /i and A, respectively: in terms of the function 

a(A) = exp(— 7riA^/4r7), 

we have, for j = 1, . . . , n, 

a(u - 2r7(/i(-''+i) + ■ ■ ■ + /i("))) ,-„Arv^-iA v" 
a{fi - 2r]{h(j) + ■ ■ ■ + h(-))) 



a{X - 2r/(/i(i) + ■ ■ ■ + h(jy)) 
The integral operator T{z, T,p, rf) is then 

T(f, r, p, 7])v = {a® Qr+p)u{z, t,t + p,r]) ® v. 



a is the operator of multiplication by the function a (A) and Qr+p is a bilinear form on 
VA[0]-valued functions, whose kernel is the Shapovalov bilinear form on Va[0]: 

®g) = j c^, ri){f{fx),g{-fx))a{fx)dfx. 

If Ai + ■ ■ ■ + A„ = 2, the universal hypergeometric function is given in terms of 
Jacobi's first theta function 6, and the phase function, see Appendix ^ by the following 
formulae: Va[0] has a basis ej = Cq ® ■ ■ ■ ® Ci ® ■ ■ ■ ® Cq, with Ci in the jth factor 
(j = 1, . . . , n) and h^^^ej = {Aj — 26ij)ej. Then the weight functions are u{t, z, A, r, rj) = 
ELi ^ait, z, A, r, r])ea, with 

. ^ _ + 2^ + ^ - - - 2^ ^.<a ^.-^ ^) n - + ^) 
' ' ' ' ^' ~ e{t - - r/A,, r) n ^(t - - r/A,, r) " 

On the other hand, uj^it, i*, /i, 77) = X] -2*! Vi v)^b, where 

v., . , _ 0{fi + 2v + t-z,- r]A, - 2r7 E^.^, A„p) " g(t - + ryA„ p) 
^''"'^'^'^^ - eit -z,- vA,,p) 9{t - z, - r^A,,p) • 

The universal hypergeometric function u is expressed in terms of these weight functions 
and the phase function 

(1 _ g27ri(2-a+jT+fcp)^^]^ _ g27ri(-z-a+(j+l)r+(fc+l)p)'j 



j,k=0 



We then have u{z, A, /i, r, p, 17) = X^a.b -^^ t", v)(^a ® e^, with 

/n 
"1 — r 
J_ _[ ^^,,A, {t - Zj, T, p)uJa{t, Z, A, r, ?7)c^f, (t, /i, J9, 77)c?t. 
, j=l 

The integral is defined as the analytic continuation from the region where Re(Aj) < 0, 
Zj e M and lm{ri) < 0. In this region, the integration cycle is just the interval [0,1]. 
After the analytic continuation, the cycle is deformed to go above the pole at Zj — rjAj 
and below the pole at zj + rjAj, as in Fig. 

The Shapovalov form is Q{X, r, ri){ea, Cb) = 5a,b<5a(A, r, rj) with 

^ g(2r/A„r)gX0,r) 

' ^' ^(A -2^ + 2^ E,<a A„ r)e{X -2r^ + 2v E,<a > r) ' 

Our first result is that our hypergeometric projective solutions of the qKZB equations 
are also projective solutions of the quantum heat equation: 

Theorem 2.1. Suppose that Ai + - ■ ■+An = 2. Let us view u{z, X, ^,T,p,ri) as a function 
u{z,T,p,r]) with values in the space of V\[0] ® V\[0]-valued functions of X and 11. Then 
u is a projective solution of the qKZB heat equation 

u{z,T,p,ri) = CT{z,T,p,ri) ® Dou{z,T + p,p,r]), 



Figure 1. The integration cycle 



where 

c = 

2TCy/Airi 

Here T{z,T,p,ri) acts on the first dynamical variable X and Dq is the operator of multi- 
plication by the function e~'^^^ of the second dynamical variable /i. 

It is then easy to construct true solutions to the system (|I]), (H) from these projective 
solutions: for any /i, l^Hj Dj{n)~^^^^Do{fiy^Pu{z, A, fx, T,p, r])C'^^P, viewed as a function 
of z, A, r, obeys (|^) and (|D in the first factor. 

In more explicit terms, we have the following statement. 

Theorem 2.2. Let for a = 1, . . . ,n 

and let the fundamental hypergeometric solution u = ^ Ua,b^a ® be defined by 

Ua,b{z, A, /i, r, p, 7]) = e~^^pb{-fx, z, p, ri)ua,b{z, A, n, r, p, rj). 

Then, for any fi,b, the V^[Q]-valued function v{z,X,t) = J2a^a.,b{z, X, fj,,T,p,ri)ea is a 
solution of the system 

v{z + p5i,T) = Ki{z,T,p,r])v{z,T), i = l,...,n, 
v{z,t) = CT{z,T,p,r])v{z,T + p), 

where C = -e^^^727rv^. 
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Proof: A consequence of (0), Theorem and the identity 

n 
i=l 

The exponential function on the right-hand side is independent of z,X,t and therefore 
does not affect the statement of the theorem. □ 

Remarks. 

1. The constant C could also be eliminated by including a factor C^/^ in u, but it is 
simpler to consider C as part of the heat equation. 

2. In the definition of the fundamental hypergeometric solution above, we have in- 
cluded an additional factor, see (||), with respect to the obvious choice. This leads 
to simpler formulae in the next section, since pb appears in the modular transfor- 
mations of the qKZB operators. 

3. Modular transformations 

The coefficients of the qKZB equations are quasi-periodic functions of zi, . . . , Zn, X with 
periods 1 and r. It is therefore natural to consider the qKZB equations as equations for 
sections of a certain vector bundle over a Cartesian power of the elliptic curves C/Z+rZ. 

Since, for any ^ ^ ^ ^ SL(2,Z), the elliptic curves with modulus r and (ar + 

b) I [cT -|- d) are isomorphic, the corresponding qKZB equations are related. For genera- 
tors, the formulae relating solutions are the following. Let for a = 1, . . . , n 

Ba[z, A,T,p,ri) = e rp ^3<k\ 3 j » Ar,T p^i^X^ z,T,ri). 

Then we have: 
Proposition 3.1. 

(i) Suppose that v{z,X) = J2a=i'^a{z, X)ea is a solution of the qKZB equations with 
parameters r + rj. Then v{z, X) is a solution of the qKZB equations with pa- 
rameters T,p,r]. 

(ii) Suppose that v{z,X) = '^^(i', A)ea is a solution of the qKZB equations with 
parameters —1/T,p/T,ri/T. Then 



n 



^(^5 A) = 2_^Ba{z,X,T,p,r])Va{z/T,X/T)ea 
a=l 

is a solution of the qKZB equations with parameters T,p,T]. 

The proof of this proposition is based on the formulae for transformation properties 
of the qKZB equations under modular transformations, see Appendix |B[ 

Now the question about monodromy is well-posed: can one express the fundamental 
solution at the transformed values of the parameters in terms of the fundamental solution 
at the original values? The answer is provided by the following result. 



Theorem 3.2. 



(i) Ua,b{z,\,^,T +l,p,7]) ^ e ^VP Ua,b{z,X,n,T,p,r]). 
(ii) Suppose that lm{r)T / p) < 0, Im(p/T) > 0. 

n „ 

Ba{z,X,T,p,r])ua,b{f,^,^,-^,f,^) =^ Ua,c{z, X, 1^, T, p, r])Mc^b{z, 1^, 1^, T, p, r])dii. 

c=l 

The monodromy matrix (M^b) is 



27ri y 477T 

The integration over /i is over the path x ^ xr] + e, x e R for any generic real e. 

The first statement of the theorem is trivial to check. The second is the first of the 
two modular relations. Introduce functions p,p^: 

Pa(X, z,T,r]) = e^'^^'^~'^^^^-^-^^^^^-^^^+^''-^'^^^^^^-'^'''^^^+^'^-^^ 

pI{X, Z,T,r]) = e^(^?+i^iA.-E?^.Ai)-^(A+27?-2r,E6 A,-2z6)(A+2r,-2r,E?+iA,)-^AE?^,A,-_ 

Then the modular relations are the identities: 
Theorem 3.3. 

(i) Suppose thatlm{r)T/p) < 0, Im(p/T) > 0. Then the universal hypergeometric function 
u satisfies the following relation, 

/ Ua,c{z,X,fX,T,p,rj)u,,b (f,-f,^,-J,-^, J 

c=l 

xQc{li,p,r))pc{-n,z,p,r))e dp, 
(6) - 2vrzy^p„(A, z, r, ^)pl f , ^) ei^ 

j<fc \j=i / 

The integration over p is over the path x xr) + e, x E M. for any generic real e. 

(ii) Suppose that lm{r)p/T) < 0, Im(r/p) > 0. Then the universal hypergeometric func- 
tion u satisfies the following relation, 

n „ 

Y j Ma,c(f, 7,7'~f'~r'r) Uc,b{z, - p, T, p, Tl) Qc{p, T, Tl) p^{p, Z, T, 7]) c'"^^' dp = 
c=l 



r/ie integration over fi is over the path x i—^ xr/ + e, x G M for any generic real e. 



This theorem is proved in 5.2 and |5.3|. 



4. The integral transformation 

Throughout this section, we assume, for definiteness, that Imr] < 0, and that rj is 
sufficiently small. The results hold for a more general range of parameters by analytic 
continuation. 

4.1. A space of functions on which the integral transform is defined. The qKZB 
heat equation is based on an integral transformation. In this section we give a space on 
which this integral transformation is defined and invertible. 

Let us fix our parameters Aj, t], t. Then the Shapovalov form has poles at the points 
—rjcTj where 

(7) a, =2(J]A,-1), j = 0,...,n, 

as well as at the translates of this points by the lattice of periods. 

At these points, the hypergeometric integrals obey the following "resonance relations" 



[V3| 



Proposition 4.1. Let r,s, 1 < a,b < n be integers and let aj = 2 Ylk<j ^fe ^ 2. Then: 

(i) If a < n, then 

Ua+l,b{z, V(^a + r + ST, fl, T, p, 7]) = e^^''^'-+'-''^+''^^+'+''^^^Ua,b{z, V^^a + r + ST, /i, T, p, 7]). 

(ii) Ui^h{z, -2r] + r + ST,fi, T,p,r]) = e2'f*'^(^i-^"+'?^i+'?^"--p)M„ 2t] + r + ST,fi, T,p,r]). 

(iii) If b < n, then 

Ua,b+i{z, X,-ri(Jb + r + sp, T,p, r]) = e^''''^"^+^-^''-^^'^+^"''^'^'>Ua,b{z, \,-ri(Jb + r + sp, T,p, r]). 

(iv) Ua^i{z, X,2i] + r + sp, T,p, r]) = e'^'"'"''^'''^~''"~'^^^~'^^"~^'^'^Ua^n{z, A, —2i] + r + sp, T,p, r]). 

Proof: (i) Using the functional relation 6{x + r + st,t) = (— 1)''"'"'^ exp(— 7ris(2x + 
st))9{x,t), we see that 

co^it, z, r/a. + r + st, t, = (^.ly^s^-^^^sit-z.^.A^y^.s^r "Q " ^ + j , 



and 

Ua^,{t, z, vcra + r + sr, r, r^) = ^.lY+s^~2.^sit-..,,-r,A.,,)-.^s^r A ^(t - z, + r^A ) 

Thus we have equahty at the level of integrands, 
(ii) Similarly, we have 

Ui{t, Z, -27] + r + ST, T, 7]) = (^_lY+s^-2^isit-z,-r^A,)-ms^r^ 

and using the relation ^ Aj = 2 we obtain 

Unit, z,2r] + r + sr, r, r]) = ^^lY+s^-^^^s(t-z^+r,A^)-.^s^r A ^(t - z, + r/A,) 

The last product of ratios of theta functions may be absorbed using the functional 
relation (^) for VL. This gives 

Un,b{z, 2?7 + r + ST, /i, T, p, 7]) = e'^^^'^'^^'^'^"^^'^ 

^ 3 

The result is then obtained by shifting the integration variable t by —p and using the 
relation 

u^,{t -p,z,f,,p,r^) = e2-(^+2,)^v(^^ 
The remaining claims (iii), (iv) are proved in a similar way. □. 

Definition. Let E^{z,T,c;ri) be the space of holomorphic functions <y9 : C — ^ ^[0], 
= Y^^aWf^a such that 

(i) If a < ra, then iPa+i{Wa + r + sr) = e'^'^'i-^^a+i-za+riAa+i+riAa) ^^^(^^^^ + r + sr). 

(ii) <^i(-2r/ + r + sr) = e2^^"(^i-^"+'?^i+'?^"+^)(^„(2r/ + r + st). 

Let E{z, T, c; rj) be the space of functions (p G E^{z, r, c; rj) such that 

(iii) There exist constants Ci,C2 > (depending on ip) such that 

|<^a(A)| <Ciexp (n^^^^ + C^lXl 



Im T 



for all a = 1, . . . , n. 



Examples of functions in E{z, r, c; 77) can be constructed using the universal hyperge- 
ometric function: 

Proposition 4.2. //Imc < then the function 

A f-^ ^ Ua,b{z, A, fi, T, -c, r])ea. 
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belongs to E{z,T,c]ri) for all values of the remaining parameters. //Imc > then the 
function 

/i Ua,b{z, A, /X, c, r, -r])eb, 

b 

belongs to E{z, r, c; 77) /or all values of the remaining parameters. 

Proof: Prop. implies that these functions belong to E^{z,T,c;ri). The bound (iii) 
follows from Lemma [4.7| below. □ 

Lemma 4.3. Lp{X) G E^{z, r, c; t]) if and only if e''*'^^/'*'''(/9(A) G t,c — T;ri) . 

Proof: It is clear that e 0(A) obeys (i) in the definition of E^{z, t,c — r^rj). Property 
(ii) follows from the identity 

for r, s G Z. □ 

For the construction of the integral transform, we will need the following variant of 
Prop. O. 



Corollary 4.4. The function 

X^e""' 4" 2^Ua,biz,\,ii,T,p,r])ea, 

a 

belongs to E^{z, t,t — p; rj) for all fi. The function 

^^e""' 4^ 2^Ua,b{z,X,^i,T,p,-r])eb, 

b 

belongs to E^{z,p, t — p;ri) for all A. 



Proof: An immediate consequence of Prop. [4.2| and Lemma |4.3| . □ 

4.2. The definition of the integral transform. We define two integral transforms 
J-'z,p,T and fiz,T,p, then show that they are inverse to each other. The first is (/? 1— *• 
^z,p,t{'-p) = <^ with 

(8) (pc{[i) = „ 2 / y(pb{-i^)Qbii^,r,r])ub^^{z,iy,ij,,T,p,-T])e'^'*vdiy. 



The second is ^z,t,pW = with 



c=l 
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These transformations depend on the choice of a real parameter a or a. We say that a 
is an admissible shift for JF if it does not he in any interval [Im(2?7 + sp), Im(— 2?] + sp)], 

s e z. 

We say that a is an admissible shift for JF if it is a generic real number. 

Theorem 4.5. Suppose that Imr] < 0. If f & E{z,t,t — p;r]) then the integral (j^ zs 
absolutely convergent, independent of the choice of admissible shift and defines a linear 
map T^^p^T : E{z, t,t — p;ri) E{z,p, r — p;ri). If ip E E{z,p, t — p^rj) then the integral 
(j^ is absolutely convergent, independent of the choice of admissible shift and defines a 
linear map ^z,t,p '■ E{z,p,T — p; rj) E{z,t,t — p; rj) . Moreover 

z,p,T ^ •F z,T,p ^^E{z,p,T—p;ri)} z,T,p ^ ^ z,p,T '^^E{z,T,T—p;ri)} 

The proof of the first part of the theorem is contained in the next subsection. The 
proof of the inversion formula is deferred to 5A . 

This theorem implies a completeness result for the qKZB equations. For generic 
^0 g (^n^ j^g^y consider the qKZB equations on the set z^ + (pZ)". Any solution is 
uniquely determined by its initial condition at A class of solutions is given by taking 
linear combinations of components of the fundamental hypergeometric solution: Let us 
say that a solution v is of hypergeometric type if it is of the form: 



Va{z,X)= / ^Uab{z,X,fi,T,p,r])Fb{fi)dfi, 



b 

for some functions Fb{fj.) and some cycle 7. 

Corollary 4.6. Any solution v{z, A) with initial condition in e ^■^ E{z, t,t — p; rf) is of 

hypergeometric type. 

More precisely, let z^ G C" be generic. Suppose ip G E{z^ , t,t — p; rj) and let (p{f^) = 
J^zo^p,r{^){n). Then, for all z e z° + (pZ)", = IlLi viewed 

as a function of fi, belongs to E{z,p, r — p-.rf) and the function v{z, X) = J2 Va{z, X)ea, 
with 

Va{z,X)= / '^Ua,b{z,X,n,T,p,r])Qb{n,p,r])e~''''^ipb{z,-n)dii, 

b 

is, for any generic 5" G M, a solution of hypergeometric type of the qKZB equations 

■hi 

with modulus t, step p and initial condition v{z^,X) = e^^^(p{X). Moreover, for all 
z E z^ + (pZ)", f (i*, A) is independent of the choice of a, and e~™^ ^^^v{z, A) belongs to 
E{z,T,T-p; T]). 

Proof: Let dj^ai.fJ') denote the eigenvalues of Dj{fi): 
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They grow at most exponentially as ^ oo. The coordinates of ijj are then 

n 

By the theorem, ip G E{z^^p,T — p;T]). Then ip clearly obeys the bound (iii) of the 
definition of E{z,p,T — p^rf). The resonance relations (i),(ii) are checked by inserting 
the definitions. For example, for a < n we have (using the fact that {z — z^)/p has 
integral coordinates) 

i)a+i{z,ri(Ta + r + sp) = ^pa{z,r]aa + r + sp)e'^''''^'°+'-^°^'^^^+'^'^^''^ 

dj^a+ii-Wa - r - sp) 



X 

3 



n 



dj,a{,-Wa - r - sp) 



On the other hand, the expression in curly brackets is equal to e^'^**^ if j = a, to 
^-2msp if j = a -\- 1, and to 1 otherwise. It then follows by Theorem [4.5| that tplz, A) = 



g-i7rA2/4»7^^^^ j)^^ is in E{z,T,T — p] T]) and that the initial condition tp{z^) = v^(A) is 
satisfied, f is a solution of hypergeometric type since the z, A-dependent part of the 
kernel of integration is 1 ® Hi -Di(/i)~^*/^M which is equal to the fundamental solution u 
up to z, A-independent factors, cf. (|^). □ 

Remarks. 

1. The result may be expressed in the following terms: our generalized Fourier trans- 
form maps the qKZB equations to the difference equations 

i>{z + p6i,fi) = Di{-fi)~^ip{z,ij,), 

which are easily solved, since Di is a diagonal multiplication operator. So the 
solution given in the corollary is the Fourier transform of the solution of this simple 
system of equations with initial condition given by the inverse Fourier transform of 
the given initial condition. 

2. An initial condition at z^ uniquely determines a solution only on the set F = 
z ° + (pZ)", since solutions defined for all zcan always be changed by multiplying by 
p-periodic functions without changing the initial condition. However the formula in 
Corollary O gives a solution for all z. On the other hand it is only for z G z^+{pZ)"' 



that we know that the integral is independent of the choice of a. For general z 
there is no cancellation of residues at pairs of poles that would allow us to move the 
integration contours, and one should expect that the solution depends on a. In fact, 
if we move a just a little bit, then the integration contour crosses infinitely many 
poles (for generic rj), so one would expect that the solution depends discontinuously 
on a. It would be interesting to understand this dependence in more detail. 
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4.3. Proof of Theorem |4.5| . We prove here that the integral transformations are well- 
defined on the considered spaces. 

We start by estimating the integrands. 

Lemma 4.7. 

I ni^ ^ M ^ ^ f (ImA)^ ^ 1^ (Imu)^ ^ 

|e ^^Uah\z, X,fi, T,p,'r])\ < Ci exp vr— h C2|ImA| + n— h C2|Im/i| 

\ Im r Im p 

for some Ci, C2 > depending on z, T,p, rj, A. 

For every e > 0,z, T,p, rj, A there exist constants C3, C4 > such that if the distance 
between A and the singularities of Qa is at least t, then 

|<9a(A,r,r/)| <C3exp f-27r2j^ + C4|Im A 



Im r 

Proof: The first bound is obtained by applying the estimates of Lemma |(J.1| to the 



integral The second follows using the lower bound in Lemma |C.1| . □ 



In order to show that the integral is independent of the choice of admissible shift and 
to bound the integral transform, we will need to shift the integration contour. The next 
result shows that this is possible since poles occur in pairs with opposite residues. 

Lemma 4.8. (i) Suppose tp = ^^a^a ^ E{z,p,T — p^rj). Then for all r,s G Z, 
c = 1, . . . ,n and a = 1, . . . ,n — 1, 



res^=-,?<7„+r+spe Uc,a{z, ><, fJ',r,p,r])Qa{n,p,v)i^a{-fi) 



-res 



■ A. -fl-i 



lJ,=-ria, 

res^=_2^+r+spe 4„ Mc,„(z, A, r,p, r7)Q„(/i,p, ?7)?/'„(-/i) 

= -res^=2r,+r+spe ''^ 4" Uc,iiz,\,iJ,T,p,r])Qi{n,p,r])ipi{-n). 

(ii) Suppose = Yl ^a^a & E{z, t,t — p;ri). Then for all r,s G Z, c = 1, . . . , n and 
a = 1, . . . ,n — 1, 

reSA=-,,a,+r+sre'' 4^ Ua,c{z,X,H,T,p,-r])Qa{\,T,r])ipa{-X) 

= -reSA=-,,<x,+r+sre''* Ua+l,c{z,\,fX,T,p,-r])Qa+l{\,T,r])Lpa+l{-X), 

resA=-27?+r+sTe'' m„,c(?, A, yU, r, p, -r])Qn{\ r, ?7)v9„(-A) 
= -resA=2,7+r+sre'' 4^ Mi^cli*, A,/i,r,p, -r7)Qi(A,p,r7)v9i(-A). 

Proof: A straightforward calculation. □ 

Proof of the first part of Theorem \4 ■ d{ . It follows from Lemma |4?^ (iii) that the integrand 
in has only simple poles at ±2r7+Z+rZ. If a is an admissible shift then the integration 
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contour stays away from the singularities and we can use the estimates of Lemma |477 . 
The integrand is then the product of a function growing at most hke e'"''^' times e^^'^ 
which for Im < converges very rapidly to zero in the real direction. So the integrand 
is an function. Moreover by Lemma (iv) the residue at 27] + r + st is opposite to 
the residue at —2ri + r + sT. It follows that the integration contour can me moved across 
each of these pairs of singularities without changing the value of the integral. This shows 
that the integral is independent of the choice of admissible shift. 

Let us next show that (p defined by (||) belongs to E{z,p,T — p^rf). Properties (i) 
(ii) follow from Corollary [4.4| , and we are left with the proof of (iii). For this we shift 

variables: let u\{z, X, fi,T,p,ri) = e'^'^Ua^biz, X, fi,T,p,ri). Then 



1 r " .^2 

'^^(Z^) = / M-t^ + l^)Qbij^ - /i, r, ri)ul^{z, u - fx, fx, t, p, -r/)e''*^c/z/. 

The integration contour is moved by this change of variable, but using the independence 
on the choice of a we may choose it to lie in the interval [0, Im r] say. 
Then Lemma 14.7] yields 



The triangle inequality — i/| < |/x| + |z/| and the Gaussian integral over u give the 
estimate (iii) in the definition of E(z,p, t — p^rf). 

The proof that @ defines a function in E{z, t, r—p; rf) is analogous with a little differ- 
ence. Now the integration contour is parallel to the line 2?7M. So if the shift a is generic it 
will not meet the singular set 5* = {±2?7} + Z-|-pZ of the integrand. However it will come 
arbitrarily close to these singularities. Still, the integral is absolutely convergent, since 
the distance to the singular set decreases polynomially: dist(/i, S") > const(l + 
for some a > 0. This implies that the divergence coming from the poles close to the 
integration contour is at most polynomial. This does not spoil the integrability which is 
due to the exponential decay of exp(— z7r/i^/4?7). □ 



5. Calculations 



5.1. The heat equation. Here we prove Theorem ^A. 
The statement of the theorem is 

„47ri»7 

Ua,b[z,X,u,T,p,ri) = —e ^ ' 



2'^^^/4M] 



(10) X 



Ua,c{z, X, fx,T,p + T, ri)Qc{fx, T + p, r])e'~ Uc,b{z, -fi, i^,T + p, p, r])dfx. 



c=l 



We proceed to evaluate the right-hand side. The /i-dependent part may be simplified 
by using the following identity: 
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Lemma 5.1. 

n 

^ z,n,T + p, r])Qc{n, T + p, r])uJc{s, z, -/i, t + p,t]) = 



c=l 



9{fi-2r] + t- s,T + p)9'{0, r + p) A g(t - Zj + r]Aj, r + p) 
e{n -2r],T + p)d{t -s,T + p) Al _ zj - r]Aj, T + p) 

e{fx + 2ri + t- s,T + p)9'i0, r + p) ^ e{s - zj + r/Aj, r + p) 
e{n + 2r],T + p)e{t - s,T + p) e{s - zj - 7]Aj, r + p)' 

Proof: Consider each term in the sum on the left as a function of fi: it is periodic with 
period one and as /i is replaced by fi+r+p, it is multiplied by exp(— 27ri(t— s)). The poles 
of the term labeled by c are at /i = /ic = 2?] — 2// J2j<c-^j at = /Xc — 2riAc = fic+i 
modulo Z+ (r + p)Z. The most general form of a meromorphic function of fi with these 
properties is 

^ ejfi- fic + t- s,T + p)e'{o,T + p) ^ ^ 9{i2 - + t-s,T + p)e'{o, r + p) 

6{fi - fic,T + p)6(t - s,T + p) ^ - ^c+i,r + p)6{t ~ s,T + p) 

The coefficients are determined by comparing the residues at the poles: 

Ac = uj^{t,z,Hc,T + p,r])uJc{s, z, -fic, r + p, 7]) 



^ jT 6{t - Zj + r]Aj, T + p) e{s - Zj + r]Aj, r + p) 
y Sit - h - V^j, r + p) J-jL eis - zj - 7]Aj, T + p)' 

Similarly one determines which turns out to be equal to — A^+i. It follows that in 
the sum over c only two terms are not canceled and we obtain our claim. □ 
The products of ratios of theta functions in the above identity may be absorbed into 
the phase functions flriA- by means of their functional relation (0). The right-hand side 
of (|10D is then 

Wt, — 7f^ I (i'~^^^'''^^^^^a{t,z,X,T,T])uj^{s,z,u,p,r]) 



271^ Ail] 

n (f^.A, (t - z, +r,r,r + p)n,,^ {s - z,,r + p, p)) XXZ'SfS'lT 
■j=i 

n 

- n (^.A. it -z,,r,r+ p)n,^^ {s -z,+p,r+ p, p)) ^v/xt;?-!!^;?^ 



dt ds d/i. 



The next thing to notice is that if we change variables in the first of these two terms by 
replacing t by t — r, s by s + p and /i by /i + 4?] we obtain exactly the second term up 
to a sign! Indeed, the shift of t in cUa produces a factor e^'^*(''*+^^) , the shift of s in uj^ 
produces a factor e~27n(jy+2»7) additional e^'^*^ from shifting the argument 
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t — s in the ratio of theta function in the square bracket. These factors are canceled by 
the shift of /i in the exponential function. 

To compute this integral we therefore have to carefully consider the deformation of 
integration contours involved in the change of variables. As we shall see, this deformation 
produces a residue at the pole t = s. 

For these considerations we assume that lm{riAj) > N, j = 1, . . . ,n, for some > 
large compared to r, p, and that the points zj are on the real axis. The general case can 
then be obtained by analytic continuation. In this range of parameters, the integration 
cycles for the t and s integration is the interval [0, 1]. The integrand, viewed as a function 
of t or s, is then regular in the strip Im(t),Im(s) G {—N,N). The first term in (|TT]), 
however, has additional poles at t = s + a + (3{t+p), (a, (3 E Z). To deal with these poles 
we move slightly the s-integration cycle into the upper half plane. After the change of 
variable, in the first term in (|ll|), t = t + r, s = s — p, the new variables are integrated 
over t G r + [0, 1], s G —p + ie+ [0, 1] for some small e > 0. The first term becomes then 
equal to the minus the second term in (0) after deforming the integration cycles to the 
original position, but during this deformation we encounter the pole at t = s. Therefore 
we obtain a residue 

(13) = yl^^"^^''/ e~^''^~''^^'^^uJa{s,z,X,T,ri)uj^{s,z,u,p,ri) 

n 

X JJ ^^.,A, {S - Zj,T,T+ p)fi^A, (S - Zj+p,T+ p, p)ds dfl. 

Using the identity (^6]) and 

— ^ / e"3^(^-^+^)'c//i = 1, Imr/<0, 

we see that this expression reduces to u{z, A, u, T,p, rj), completing the proof of Theorem 



5.2. First modular relation. Here we present the proof of Theorem p.3| (i). The 
assumption that Im(?7r/p) < implies that the integration over /i on the path x ^— T]x+e, 
X G R, converges absolutely: the singularities of Q are avoided for generic e, and at 
infinity the factor exp{—iTTT^'^/4rip) converges very fast to zero. To get the identity as 
stated in the Theorem, one uses the following simple properties of pa- 

rfe\^p'p' p' p) \ p' p' p' p J 

Let US proceed with the proof. We insert in the left-hand side the hypergeometric 
integral for ^ and Uc^h and call the integration variables t and s, respectively. It will 
be convenient to make the change of variables s ^ s/p. To apply Lemma we use the 
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transformation properties of weight functions: 

A similar formula involving gives the transformation behavior of . Using Lemma 



bA\ and the functional relation (|3l|) of fia, we get: 

n „ 



? /£ £ 1 T 2 

) ' p'' p'' p ' P ' P 



C=l 

(13) xQ^{fi,p,r])pc{-fi,z,p,r])e~^^ 



X 

where 



- / uJa{t,z,\,T,r])uj'^{s/p,z/p,u/p,-T/p,r]/p)e ^"^^ '"/^^^ 



' ^ " ^(/i-2r^,p)^(t-.,p) 



As in the proof of the heat equation, the two terms in this equation cancel after formally 
changing variables th^t — t, s^s — t, fih^fi + At] in the first term. However, 
we have to carefully see what happens to the integration cycles after this change of 
variables. We consider the region of parameters where zj G M, Im(?7Aj) >> and 
7r/2 >> arg(p) > arg(r) > 0. 

Then the original integration contours may be chosen to be the interval [0,1]. For 
z E C, let us denote by 7^ the path 

r rz, r G [0, 1]. 

After the change of variable s —>■ s/p, the s integration contour becomes the path 7^. 

Lemma 5.2. Suppose that p, r are complex numbers in the upper half plane. Let f{t, s) 
be a meromorphic function of two variables such that f{t + 1, s) = /(t, s + p) = f(t, s) 
and such that a{t, s) = {t — s)f{t, s) is regular in a domain containing 

{{t + rr, sp + rr) 1 1, r, s G [0, 1]}. 

Then 

I f(t + T,s + T)dtAds= / f(t,s)dt Ads + 2TTi / a{t,t)dt, 
and all integrals are absolutely convergent. 
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Proof: Since p is not real the integration contours intersect transversally, so that f{t, s) 
(which has at most a simple pole on the diagonal) is absolutely integrable. More gener- 
ally, with our assumption on the regularity of /, the integral 



f{t + rr, s + rT)dt A ds 

'7iX7p 

is absolutely convergent for all r G [0, 1]. Let, for e > 0, be the integration domain 
obtained from 71 x 7^ by removing the points where |t — s| < e and set 



-^r-(e) = / f(t + rT,s + rT)dtAds. 



Then Ii — Iq = lim^^o /q j;:^ri^)dr. We have 
d 



/^(e) = -r / d{f{t + rT,s + rT){dt- ds)) 

dr 



—T I f{t + rT,s + rT){dt — ds) 



— r 



dt — ds 

a[t + rr, s + rr) 

9_D< t — s 



= T a{u + rT,rT) \- 0{e) 

J\u\=e U 

= 2iTiTa{rT,rT), 

from which the claim follows immediately. In the application of Stokes' theorem above 
the other boundary components of dD^ do not contribute since the function is periodic. 
The remaining component is homotopic to a circle with negative orientation, which 
explains the change of sign. □ 
Applying this lemma to our situation after shifting fi by Arj in the first term yields: 
f 7700 



(HD = - / / ^a{t, z, A, r, T])u^{t/p, z/p, u/p, -t/p, r]/p)e'^^^^'"^^'^'''^ 

n^^.A,(t - z„p,r)n.^^ -1, rfte-^-^rf/z. 



Now the Gaussian integral over /i along r)M. may be performed explicitly, and our claim 
follows from (^) and 

5.3. Second modular relation. Here we present the proof of Theorem |3.3| (ii) which 
is parallel to the proof of Theorem |3.3| (i). 

The assumption that Im(?7p/r) < implies that the integration over fi over the path 
X ^— + e, X G M, converges absolutely: the singularities of Q are avoided for generic 
e, and at infinity the factor exp{—iTTpfi'^ /Arjr) converges very fast to zero. 
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We insert in the left-hand side the hypergeometric integral for Ua,c and Uc,b and call the 
integration variables t and s, respectively. It will be convenient to make the change of 
variables t —>■ t/r. To apply Lemma |5.1| we use the transformation properties of weight 
functions: 

(14) f, f , = pr(/^,^^,r,ry)-^ev(2-E.-.A.)(/^+2^).;V(t,£,/.,r,ry). 

We have 



IL n 

i\ \/ i'K'p 2 

c=l 

= ^ j UJa{t/T,z/T,X/T,-p/T,ri/T)uJ^{s,Z,-U,p,Vl)e~'^^^^^''''^^ 



where 



<f>(nts) - ^(/^ - 2r/ + t - r)g^(0, r) z.(2t-y:;-_, .,A,)(;.+2r,) 
^ " e{p-2rj,T)e{t-s,r) 



X 

3=1 

The two terms in the right hand side of this equation cancel after formally changing 
variables t^^t — p, s^s — p, fi^^fi — At] in the second term. However, we have to 
carefully see what happens to the (t, s) integration cycle after this change of variables. 

We consider the region of parameters where Zj G M, lm{riAj) >> and n/2 » 
arg(r) > arg(p) > 0. In this case the {t, s) integration cycle is the product of paths 
7t X 7i. 

Applying to this situation Lemma |5.2| after decreasing p by Arj we see that the left 
hand side of the new modular relation is equal to 



27ii 

T 



rjoo 



rjco J 7p 



UJa{t/T,z/T,\/T,-p/T,r]/T)uj'^{t,Z,U,P,T])e 2" ^^^^ "''^6 4^^'''' 



j=i \ / 

Now the Gaussian integral over /i along rjM. may be performed explicitly, and our claim 
follows from (^) and i\L2^ . 
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5.4. The inversion formulae. Here we conclude the proof of Theorem [4.5| . We have 
to prove the identities 



(15) 



^ Uc,aiz, -/i, A, p, T, -r])u^^b{z, ij,, ly, p, r, r])Q^{~iJ,, p, r]) dfi 



SabS{X + v) 



(16) 



1_ f ^ 

We give a proof of the first of these identities. The second is treated in a similar way. 

As in the proof of the heat equation, the /^-dependent part of the integrand may be 
simplified by collecting the terms with the same poles, and using the functional relation 
of i^a, with the result: 

1 /■ " 

n / y]«c,a(2', -/i, A,p,r, -77)uc,fe(z,/i, z/,p,r,?7)Qc(-/i,P,??)c?Ai 

c=l 

ujlit, A, r, -r/)c.,^(., r, r/)e-^(^+^)^ 

g(/i + 2r7 + g-t,p)g^(0,p) fr^ 
(17) - ^(;, + 2r,,p)^(.-t,p) - %^P^-)^.A. (^ - 

g(/i-277 + g-t,p)gXO,p) -,^n ^ r ^ ^ 



(f. A, /i, r, p, r\)UbXz, v, -/x, r, p, -r])Qc{lJ,, p, v) djj, 



^ e{fi- 27], p)e{s-t,p) 

As in the proof of the heat equation, the two terms in this equation cancel after formally 
changing variables t \—>- t — r, s ^ s — t, fi fi+Arj in the second term. However, we have 
to carefully see what happens to the integration cycles after this change of variables. 
We consider a region of parameters where Zj G M and rj, r are small compared to 1, 
Im(p) and to the distance between the zj. We study the singularities of the integrand 
and the form of the integration cycles in the vicinity of zj. The poles of the integrand of 
u{z, /i, z/, p, T, rj) are at s = Zj — rjAj — kr — lp and at s = Zj + rjkj + kr + lp {k, I G Z>o). 
The cycle for the integration over s goes below Zj + r/Aj and above Zj — fjAj. The poles 
of the integrand of u{z, —fi, A, p, r, —rj) in the vicinity of Zj are at t = zj + rjAj — kr — lp 
and at t = Zj — rjAj + kr + Ip {k, I G Z>o). The t integration cycle goes above Zj + rjAj 
and below Zj — rjAj. These integration contours are depicted in Fig. ^. 

To treat the two terms on the right-hand side of (^^ separately, we split the contour 
for the t integration into two pieces, as shown in Fig. ^. Then the integration cycle does 
not meet the additional poles at t = s of both terms. 
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+ 

X 




+z,-riA^-x 

X 

Figure 2. The integration cycle. The sohd hne is the contour for the 
integration over t and the dashed hne is the contour of integration for s. 
They are oriented from left to right. The symbols x, + indicate the poles 
of the integrand as a function of t and as a function of s, respectively. The 
points Zj ± r/Aj are poles for both variables 

+ 

X 

'(Sj-nAj""" 

X 

+ z,-riA^-x 

X 

Figure 3. The integration cycle of Fig. |^ can be deformed to this cycle. 
The circle around zj — rjAj is oriented counterclockwise 



Let us first consider the first term in (0). As a function of t it is actually regular 
at Zj — rjAj so that the circle around this point does not contribute. In Fig. the 
remaining integration cycle is shown, along with the position of the poles for t and s. 

After changing variables t \—>- t — t, s^s — t, n^fi + Ar] in the second term of 
(p!7|) the integrand of the second term becomes equal to the first one, but the integration 
contour are shifted: the /^-integration cycle can be deformed back to the original one 
without encountering singularities. The cycle for the t and s integration is depicted in 
Fig. |. 

We have to deform the integration cycle to the position shown in Fig. ^. In so doing, 
we pick the residue at t = zj — rjAj + r of the residue at s = t. At this point, we may 
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X 



X 



X 



Figure 4. In the first term of tlie riglit-fiand side of ([17D the integrand 
is more regular and the cycle may be replaced by this one 



+ 

rx 



X 



X 



Figure 5. This is the integration cycle after the shift of the variables t 
and s in the second term. The poles of the result are indicated by the 
symbols x , + as in the preceding figures 



use the identity (^Ol). The result is 
-(27rz^2 



-rest= 



reSi=2^._^A,t^a ^) ^) -V)^bi't, z, u, r, r])S{\ + v). 



The residue may be computed at zj — rjAj since its argument is r-periodic. Our claim 
follows then from the 

Lemma 5.3. 

n 

^ rest=^^._^A.u;^(t, z, A, r, -r])uj^{t, z, -A, r, r]) = 5abQa{.-^, r, r])'^. 



Proof: This formula can be deduced from the more general result in [[rV |, formula 



(C.4). For the sake of completeness, we include the simple proof in this special case. 
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The product of weight functions of which we compute the residue is an elhptic function 
of t whose simple poles in a fundamental domain are in the set {zj ± riAj,j = 1, . . . , n}. 

If a > b, then the function is regular at zj — rjAj, j = 1, . . . ,n, so the residues vanish. 
If a < 6, then the function is regular at Zj + rjAj. Therefore the sum of the residues in 
the claim is the sum over the residues at all poles and vanishes by the residue theorem. 

If a = 6, only the residue at t = Za — rjAa gives a non-vanishing contribution, which is 
easily evaluated, and gives Qa{—X,T,ri)~^. □. 



6. SL(3, Z)-IDENTITIES FOR HYPERGEOMETRIC INTEGRALS 

In this section we recast our results into a form which shows the analogy with the 



identities discovered in FV2 | for the elliptic gamma function. In |FV2] we showed that 
the elliptic gamma function is a "degree 1" generalized Jacobi modular function for 
the group G = SL(3,Z) k Z^. This amounts to the identities of Theorem below. 



The hypergeometric integral u obeys a non-Abelian version of these identities. We also 
discuss the geometric interpretation of these identities: they can be interpreted as the 
projective flatness of a discrete connection on a vector bundle over a suitable G-space, 



see also ||FV4|] . As in the case of Gamma functions these properties are more transparent 



if we rewrite the identities, which we formulate in in homogeneous variables x defined 
by T = X1/X3, p = X2/X3. The G-action is then linear. The next step, as in the case 
of r, is to extend the range of parameter to a G-space, as the condition of positivity 
of imaginary parts is not preserved by the action. This can be done in a surprisingly 
easy way by reflection arguments. We do it in the simplest case of one tensor factor 
in |6.5| . After this is done, we rewrite the identities as the projective flatness of an 
SL(3, Z)-connection in B^. 



6.1. Identities for the elliptic gamma function. We first recall the properties of the 
elliptic gamma function that are relevant here. More details are included in Appendix 

E 

Theorem 6.1. |p^'Vl|| The elliptic gamma function, defined by the formula 

^ 1 _ p27ri({j + l)T+(fc + l)p-2) 

r{z,T,p)= n ^ 

j,k=Q 



o2-Ki(jT+kp+z) 



/or Im(r), Im(p) > 0, obeys the following identities 



(18) 
(19) 
(20) 



r(z + l,r,p) = T{z,T,p), 

T{z + p,T,p) = eo{z,T)T{z,T,p), 

T{z + T,T,p) = eo{z,p)T{z,T,p), 
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(21) r(z,r + l,p) = T{z,T,p + l) = T{z,T,p), 

(22) Tiz,T + p,p) = ^("'^-P) 

r(z + T,T,p + T) 

(23) r(^/r,-l/r,p/r) = e^'^«(^^^'^)r((2 - r)/p, -r/p, -l/p)r(;., r,p), 

for some Q G Q(r, p)[2;] given in Appendix^. 

In fact r is also defined for negative imaginary parts of r and p by a reflection proce- 
dure, see |A.5|, and the above identities continue to hold in this wider range of parameters. 



6.2. Identities for hypergeometric integrals. Let Es{z, r, c; rj) be the space of func- 
tions Lp G E^{z, T, c; rj) such that there exist constants Ci, C2 > (depending on (p) such 
that 

I '^'^^^ /AM ^ ^ f (ImA)^ ^ , , 

|e </^a(A)| < Ciexp 7ri- ^ + C2IA 

y Imr 

for all a = 1, . . . , n. 

For 5 = we have constructed examples of functions in Es in Prop. |4.2|. 



Proposition 6.2. Lei ^^(z, r, c, r^), $t(^) i/ie operators acting on V/<^[0]-valued func- 
tions ip{X) = ^^^^ (PaW^a as 

{^s{z,T,C,ri)ip)ai\) = e Pa{X,Z,T,r|)(pa{- 



T 



Then these operators restrict to isomorphisms 

(Z \ C 71 \ 

-, ^ E6+c{z,T,C,T]), 
T T T T J 

^t{v) ■ Esiz, r, c, r]) Es+i{z, T,c+T,r]), 

■ Esi-z,T,-c,-r]) ^ E_siz,T, 0,7]). 

Then Theorem about the Fourier transform may be reformulated, in a slightly 



generalized form, as follows: 
Proposition 6.3. The operator 

/n 
^ Uab{z, A, /i, r, p, r])Qb{fi, p, r])ipb{-fj.)dfj., 
b=i 

is an invertible linear map from Es{z,p, r, rj) onto E_i/s{z, r, —p, rj). 
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The properties of the universal hypergeometric function may then be expressed as 
relations for these operators: first of all the analogue of the identities (p!8| ) -([20| ) are 
the qKZB and mirror qKZB equations (^. The remaining identities are given by the 
following result. 

Theorem 6.4. (i) Es{z,t + r,p + s,ri) = Es{z,T,p,ri) , if r, s E Z,. The identities 

V{z,T + l,p,r]) = V{z,T,p + l,r]) = V{z,T,p,r]) 

hold m }lom{Es{z,p,-T,T]),E_i/s{z,T,p,r])). 
(ii) The identity 

V{z, T, p, 7]) = Ct ^t{v) 'r,r + p, r]) $r(?7) V{z, T+p,p, r]) ^t{v) 
holds in B.om{Es{z,p,r,ri), E_i/s{z,t, —p,!])) with 



Ct 



271^ Air] 
(iii) The identity 

^s{z, r, -p, r]) V (I, -i, f , ^) $5 (I, £, -i, ^) 
= cs V{z, T,p, ri) ^siz,p, T, ri) V (^|, -i, -j, 2) 
holds in B.om{Es{z/p, —r/p, —l/p,ri/p), Es'{z, r, —p, r])) with 5' = p5 /{I — t6) and 



Proof: The first statement (i) is trivial. The second statement (ii) is a reformulation of 



Theorem 2.1 and fiii) is a reformulation of Theorem 3.3. □ 



6.3. The universal hypergeometric function for one tensor factor. For simplic- 
ity we assume from now on that the tensor product Va of 5^2 representations consists of 
one factor which is the three-dimensional representation V2. In this case the universal 
hypergeometric function does not depend on zi and is a scalar function, 



iirXfi 



e{x + t,T)e{fi + t,p) 



(24) u{K^.,r,p,,) = e- / ^^.(t, r,p) ^^^ _ ^^'^ ^^^^^ _ ^;^^^^ cit. 



c 



The integration cycle C is as in Sect. |^: in this case it can be defined to be the interval 
[0, 1] if T] has positive imaginary part. For general 77 the integral is defined by analytic 
continuation, with the effect of deforming the contour as in Fig. |l|. The resulting function 
is a meromorphic function on C x C x x i7+ x C, where denotes the open upper 
half-plane. This function is regular as long as the integration cycle is not pinched between 
poles of the integrand. As these poles are at ±(2?7 + rr + sp + t) , r,s G Z>o, t G Z, m is 
regular if 4t] does not belong to the set Zr + Zp + Z. 
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The following is an interesting alternative formula for the function u, 



X J2 e-2'^*(^'^+'^^)+(2^'+^)(2'=+^)''^o(A + 2ri + kp, r)eo(/i + 2r/ + jr,p) 

j,k=0 

9o{lp + 4r],T) 9oilT + Ari,p) 
fj^^o((/ + l)p,r)fj0o((/ + l)r,p)' 

where 9o{z,t) = YlT=o('^ " e'^ni{jr+z)^(^i _ ^2nii{j+i)r~z)^^ ^pp^ ^ rpj^ig formula can 
be proved by moving C to infinity and picking up the residues at the poles, and is valid 
in a certain region of parameter space. We will not use this formula in this paper. 
The Shapovalov form is 

^^^'^'^^ 9{\-2r^,T)9iX + 2r^,T)- 

6.4. Spaces and operators. Consider and the projectivization of the dual space, 
P{C^)*. Consider W C {C - Of x P{C^)* where 

W = {{{xi,X2, X3), {yi : y2 : Z/s)) I a^iZ/i + X2y2 + x^Vz = } . 

The natural projection 14^ — > — is a projective line bundle. The group 5^(3, Z) 
acts g: {x,y) {gx, {g^Y^y) for g G SL{?,,Z). 

For (x; y) eW such that Ivcixi/x^ 7^ 0, ?/2 7^ introduce the space of functions 

{v{\) I v{x,X) e Es=.y,iyAf^, f|, 4)}, if Im > 0, 
(^;(A) I v{x,\) e Es=-y,/yA-f^. -t, 2^)}, if Im < 0. 

Thus E{x] y) consists of entire holomorphic functions f (A) obeying the resonance relation 
f (— 1 + rx3 + sxi) 



E{x;y) 



the bound 



(25) 



^2ms{2(T+x2)/xsyl^'^ _l_ _|_ gxi), where a is the sign of Iraxi/x^, and 
niyiX^ 



v{X)e 2x31/2 



(Im^)2 
< Ciexp I n- ^ + C2IAI 



Im^i 

X3 



for some constants Ci, C2 > 0. Note that 

E{xi,X2, X3; yi, y2, 1/3) = E{-xi,X2, -x^; -1/1,1/2, -ys). 

For Im^ > 0,lm^ > 0, introduce an integral operator U{xi,X2,X3) acting on func- 
tions in E{x2, Xl, X3; 1/2, yi, ys) as 
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The integration path is given as follows. Suppose first that \x2\, jxsl >> 1. Then 
the integration is over a straight line which does not intersect the segments joining 
— l+rx2 + sxs to l+ra;2 + sa;3 (r, s e Z) and such that the integrand decays exponentially 
at infinity along it. Such a path exists as the integrand behaves as exp(— const /i^) at 
infinity, which decays for n in a. cone. The resonance relation implies that the integrand 
has opposite residue at the pairs ±1 + rx2 + sx^, so that, up to sign, the integral is 
independent of the choice of path. For X2,X3 general, the integral is defined by analytic 
continuation. The operator is defined up to multiplication by ±1 since one needs to 
choose a square root and an orientation of the path. 
We assume that 

2^Zxi + ZX2 + Zxg. 



This ensures that f/(xi, X2, Xa) is well-defined. Indeed, the integration kernel (2^) defin- 
ing U is regular at points obeying this condition by the discussion of |6.3| above. Moreover 
this condition implies that in the integral over /i defining U, the integration contour is 
not pinched between poles. 

Proposition |6.3| says that f/(xi, X2, Xs) is an invertible linear map 

U{xi,X2,X3) : E{x2,xi,X3;y2,yi,y3) ^ E{xi, -X2,X3;yi, -^2,1/3), 
if Im 2i > 0, Im 2a > 0, 

X3 '2:3 ' 

Introduce operators a^x^), X2, xs) where 



ITT /' 1 \ 2 1\ 



I r, ( ^ \^ I \^ ill 27ri ™ 

(3{xuX2,xM^) = e^^^A^^^ r«™3^(^) iflm^>0, 

I T1 r ^ .\2 ^\ .\2 I 1 I I 2jri ™ 

P{xi,X2,X3)v{X) = e ^^^^^ iflm^<0. 

The A-independent terms in a and /3 were added to simplify the relations (l)-(4) below, 
which would otherwise only hold up to factors depending on x. 

Proposition 6.5. The operators a and (3 are isomorphisms 

a{x3) : E{x2, Xi, Xs; 2/2, Z/s) ^ E{x2, Xi + X2, X3; y2 - 2/1,1/1, 2/3), 

(3{xi, X2, X3) : E{x3, xi, -X2; y3, yi, -2/2) E{x2, xi, X3, 2/2, 2/1, 1/3), 
for all {x, y) & W for which they are defined. 

Proof: The statement for a is easily checked inserting the definitions. The fact that f3 
respects the resonance relation is also easily checked. We are left to prove that v obeys 
the bound (^) if and only if jSv does. Let us prove the only if part in the case where 
X2/X3 has positive imaginary part. The other cases are proved in the same way. So 
we assume that | exp(7r2?/3A^/2a;2yi)f (A)| < Ci exp(7rj|^^^^^^ + C2IAI). It follows using 



28 



the relation ^ Xii/i = that v' = P{xi, X2, a^s)^ obeys the bound | exp(— 7ri?/2A^/2a;3?/i + 
TciX"^ / X2Xs)v' {X)\ < C[ exp(7rj[^^^^^ + C2|A|). The claim then follows from the identity 



x2 flm— ^ flm— ^ 

Im 



X2X3 Im — Im — — 



□ 



Theorems 12.21, 13.31, 14.51 take the form of identities: 



1. The q-heat equation, 

a{x3)U{Xi,Xi + X2,X3) a{x3)U{Xi + X2,X2,Xs) a{x3) = U{Xi,X2,X3) , 

holds on E{x2, Xi, X3; ?/2, Vi, 2/3), if Imf^ > 0, Imf| > 0. 

2. The first modular equation, 

U{xi,X2,X3) P{xi,X2,X3)U{-X3,-Xi,X2) = 

(3{-X2, Xi,X3) U {-X3, X2, Xi) I3{-X3, X2, Xi) , 

holds on E{-xi, -X3, X2] -yi, -1/3, Z/2) if Imf^ > 0, Imf| > 0, Imf^ > 0. 

3. The second modular equation, 

U{-X2,-X3,Xi) P{-X2,-X3,Xi)U{xi,X2,X3) = 

P{X3, X2, -Xi) U (Xi, -X3, X2) P{Xi, X3, -X2) , 

holds on E{x2,Xi,X3;y2,yi,y3) «/Imf^ > 0,lmf| > 0,lmf^ > 0. 

4. The inversion relation, 

U {Xi, X2, X3) U {-X2, -Xi, -X3) = 1 , 

holds on E{-xi, -X2, -X3, -yi, -y2, -y3) if Imf^ > 0, Im^ > 0. 

Notice that each of the operators U in these identities is defined up to multiplication by 
±1, so the right hand side of each of the identities is equal to the left hand side up to 
multiplication by ±1. 



6.5. New range of parameters. Extend the definition of the operator U{xi,X2,X3) 
from the 
formulas 



from the domain Im^ > 0,lm^ > to the domain 0:1/0:3, 0:2 /X3 G C — M by the 



U(xi,X2,X3) = U(x2,—xi,X3) ^ if Im— < 0, Im— > 0, 

3^3 X3 

U(xi,X2,X3) = U(—X2,xi,X3y^ if Im— > 0, Im— < 0, 

2:3 3:3 

U{xi,X2,X3) = [/(— Xi, — X2, X3) if Im— <0,lm— <0. 



X3 X3 



Theorem 6.6. 
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(i) For any {x;y) G W such that Xi/x^, X2/X3 G C — M and 2 ^ Zxi + 'IjX2 + 
Zxs, the operator U{x) defines an invertible linear map i?(x2, xi, X3; 1/2, l/i, 2/3) —>■ 
E{xi,-X2,X3;yi,-y2,y3)- 

(ii) Moreover, the q-heat equation, first and second modular equations, the inversion 
equation hold for Xi,X2,X3 such that Xi/x2,Xi/x3,X2/xs G C — M and 2 ^ Zxi + 
7^X2 + ZX3. 



The six-term relations of Theorem |6.6| are the commutativity of the following diagrams. 
The q-heat hexagon: 

E{x2,xi,X3;y2,yi,y3) 

a{x3.) / \ f/(xi,X2,X3) 

E{X2, X1+X2, X3; 2/2-2/1, 2/1, 2/3) E{xx,^2, 2/1, -i/2, 2/3) 

f/(Xi+X2,X2,X3)| {^(a^s) 

E(xi+X2,-<r2, xs; 2/1, 2/1 -2/2, 2/3) ^(xi,-a;i -X2, X3; 2/1-2/2,-2/2, 2/3) 

a(x3) \ / f/(Xi,Xi+X2,X3) 

E(X1 + X2, X2, X3; 2/2, 2/1-1/2, 2/3) 

The first modular hexagon: 

E{-xi, -X3, X2; -2/1, -2/3, 2/2) 
/3(-X3,a;2,xi) X \U{-X3,-xi,X2) 
E{x2, -X3, xi, 2/2, -2/3, 2/1) ^(-a;3, a;i, X2; -2/3, 2/i, 1/2) 

[/(-X3,X2,Xi)| |/?(xi,X2,X3) 

E{-X3, -X2, xi, -2/3, -2/2, 2/1) ^(2^2, X3; 2/2, 2/i, Vs) 

P{-X2,Xi,X3) \ / U{Xi,X2,X3) 

E{xi,-X2,X3;yi,-y2,y3) 
The second modular hexagon: 

E(x2,xi, 0:3; 2/2, 2/1,1/3) 

/3(xi,X3, -X2) ^ \U{xi,X2,X3) 

E{-X3, Xi,X2] -2/3, 2/1, 1/2) ^^(2;i, -X2, X3; 2/1, -2/2, 2/3) 

f/(xi, -X3,X2)| |/?(-X2,X3,Xi) 

X3, X2; 2/1, 2/3, 2/2) E{-X3, -X2, xi; -2/3, -2/2, 2/i) 

/3(x3,a;2, -xi) \ ^ f/(-a;2, -a;3,xi) 

E(-X2,X3,Xi; -2/2,2/3,2/1) 



The proof of Theorem 6^ is done by reduction to the case where the imaginary parts 
of X1/X3 and X2/X3 are positive. In fact it is straightforward to see that in all cases the 
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identities for general imaginary parts can be rewritten using the definitions as identities 
at some other values of x, y where the imaginary parts are positive. 

6.6. Projectively fiat SL(3,Z) connections. We give here an interpretation of our 
relations in terms of projectively flat discrete SL(3, Z)-connections. 

We start by introducing the notion of discrete connections. Let X be a G-space 
with a fixed presentation by generators e^, j = l,...,k and relations Ri = 1, i = 
1, . . . ,r. A discrete G- connection on a complex vector bundle n : F —>■ X, with fibers 
F{x) = 7T~^{x), X G X, assigns to each generator cj a collection of linear isomorphisms 
(j)ej{x) : F{e~^x) — F(x), x G X. The parallel translation along an element w of the 
free group Freefc generated by the e/s is the collection of maps 4>w{x) ■ F{w~^x) F{x) 
uniquely defined by the properties 0i(x) = Idpix), (pww'{x) = 4>w{x) o (p^i^w^^x). Here 
w ^— w is the canonical projection FreCfc — > G. The curvature of a discrete G-connection 
is the collection of parallel translations 0_Rj(x) G End(F(x)) along the relations. A 
connection is called projectively flat if (pji- (x) G CId for all i and x. 

Let TT : F — s> y be a vector bundle over a subset y of a G-space X. Then a discrete 
connection defined on Y assigns to each generator cj a collection of linear isomorphisms 
(j)ej{x) : F(eJ^x) — > F{x), for all those x E Y such that e~^x G Y. The parallel 
translation 0^(x) is then defined on some subset Y^ gY. A discrete connection defined 
on Y is projectively fiat if 4>rXx) G C Id for alH = 1, . . . , r and x G Yr.. 

In our case, G = SL(3, Z) and for X we take certain orbits in W. 

The group SL(3,Z) is generated by the elementary matrices Cij, {I < i, j < S,i j) 
with ones in the diagonal and at {i,j) and zeros everywhere else. The relations can be 



chosen [Ml to be 



(26) CijCjk = eikCjkeij, z,j, distinct, 

(dsej/eis)^ = 1- 

An SL(3, Z)-orbit X in is called regular if for all (x, y) E X 

1. Xi/xs G C — M and X2/X3 G C — M and 

2. 2 ^ Zxi + Zx2 + Zx3. 

Let X be a regular orbit in W. For (x, y) G Y^ = {{x,y) G X | yi 7^ 0}, set 

y) = E{x2, xi, X3; ?/2, yi, y-i)- 
For (x, G = {(x, G X I t/a 7^ 0}, set 

F'^'ix; y) = E{xi, -X2, X3; 2/1, -?/2, ys). 



Proposition 6.7. Let X G W be a regular orbit. 
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(i) The assignment Cij (pijix) with 

(pi2{x) = a{xsy^U{xi - X2,xi,X3y^a{x3y^, 
013 (a;) = 023 (a;) = 1, 
02i(a;) = a{x3)'^, 

031 (a;) = P{-X2,Xi - X3,X3), 

032(a;) = P{X3 - X2,X3,-Xiy^U{x3 - X2,X3,-Xi)~^P{-X3,Xi,X3 - X2y^, 

defines a projectively flat SL(3, Ij)- connection defined on G X . 

(ii) The assignment Cij i— > 4'ij{x) with 

ip^2{x) = a{x3), 

i'nix) = ip23ix) = 1, 

tp2i{x) = a{x3)U{x2,X2- Xi,X3)a{x3), 

i'3lix) = -X3,-X3,X2)"^f/(Xi -X3,-X3,X2)"^/?(X3,X2,X3 

^^32(3;) = /5(a;i,X2 - X3,X3), 

defines a projectively flat SL(3, X)- connection defined on C X . 

(iii) For all i ^ j and all x eY'*' nY"/" such that t^^x ^Y'^ ^Y^ , we have 

U{x)i)ij{x) = (f)ij{x)U{e^^x). 

Proof: It first follows from Theorem ^]6| (i) and Prop. |6.5| that (f)ij{x) is in all well- 
defined isomorphism from F'^(e^-^(a;, y)) to F'^{x, y) and similarly for if). The other claims 
are then simple consequences of the q-heat, modular and inversion relations (Theorem 
|6l6| (ii)), and the relations 

U (Xi, X2 + X3, X3) = U{xi+ X3, X2, X3) = U {xi,X2, X3), 

which follow from the fact that u{X, iJ,,T,p,ri) is 1-periodic in r and p. 

The easiest way to do the computations is to first check (iii), which can be easily 
deduced from our three term relations. This identity implies that if the curvature as- 
sociated to a relation is scalar for one of the connections (p or tp, then it is scalar (and 
equal) also for the other connection. Then the curvature can be computed using or 
ip, whichever is simpler. For example, to compute the curvature Cf| = (pR associated to 
the relation R = 6126326^2^^2 ^ better to use the connection ip. One gets 

^i2(a;)^32(er2^x) = Cf2(a;)V'32(x)V^i2(e32^x), 

where 

C'(i{x) =exp 



,32 , 2m X2 



9XiX3(xi - X2){X2 - X3) ' 

□ 

In particular, if an orbit does not contain any point with ?/i = or ?/2 = 0, which is true 
for generic orbits, we have Y'f' = Y^ = X and the connections are defined everywhere. 
The trouble is that we do not know if the spaces F'^, F^ are nontrivial for these orbits. By 
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contrast, we have infinitely many examples of linearly independent functions in F'^{x, y) 
with ?/2 = or F'^{x, y) with yi = 0. Indeed, these spaces are isomorphic to spaces 
Es{t, —p,ri) with 6 = 0, which contain the functions : A i— > u{\, iJ,,T,p,ri) with 
any fixed /i, see Prop. [4.2| . These functions are non-zero for generic /j, since they are 
the meromorphic kernel of an invertible integral operator. There are infinitely many 
linearly independent functions among them since they behave differently under shifts 
by 1: u^{\ + 1) = -e-^^^/2ry^^(A). In this case we may construct a connection defined 
everywhere by gluing the two partially defined connections: 

Theorem 6.8. Let X be a regular orbit containing a point {x,y) with yi = 0. Let Y^, 

G X be the subsets as above, on which the connections (p, ip are defined. We have 
Y't'ijY'^ = X. Let F be the vector bundle on X obtained from F^ and F^ by identifying 
the fibers over 

yi' n via U{x) : F^{x,y) F^{x,y). Let a connection Xiji-^ 

) be 

defined as (pijix), if both {x,y) and e~^{x,y) are in Y'^, and as tpijix) if both {x,y) and 
e~j^{x,y) are in Y"^ . Then x a well-defined projectively flat connection on X. 

Proof: Since any two points {x, y), e~j^{x, y) E X related by a generator belong both to 
Y^ or both to Y"^ , the connection Xij is defined in all cases. Moreover, by Prop. (iii), 
X is well-defined on F, at least up to sign. It is easy to check that that the curvatures 
along the relations involve products of isomorphisms Xij{^) niapping between fibers at 
points {x,y) which are all in Y'^ or all in Y^ . Thus the claim that x is a projectively 
fiat connection follows from the projective fiatness of (f) and ip. □ 

6.7. Formula for the curvature. We give here a formula for the curvature of the 
discrete connection (and ip). It is convenient to do a gauge transformation (pij^x) = 
(pij{x)gij{x), with 

^12(0;) = exp ( — ^ — ), ^13(3;) = exp ( — ^ — ), 

\9 xi [xi - X2) X3 J \9xi (xi - X3) X2/ 

^31 (x) = 1, ^32(x) = exp ( - ) . 

V 9 X3 {X3 - X2) Xi J 

Then the curvature of the connection has components C^-{x), C{x) associated to the 
relations (^61) . By definition, they are given by 

(t)^j{x)(j)kl{e-^x) = Ci-{x)(f)ki{x)(f)ij{e^lx), i^l, j ^ k, 
(j)ij{x)(j)jk{eij^x) = Cl^{x)(f)ik{x)(f)jk{eikX)(f)ij{eJ,^e^,lx), ij, k distinct, 
4>s{x)(l)sis~^x)4>s{s~^x)(l)sis~^x) = C{x), 
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where sx = {—X3,X2,xi) and 

(psix) = (i)i3{x)(j)3i{e3ieiixy^(j)i3{e^^e3ix). 



We have 



Cr^ix) = Cllix)-' = exp 



27li X2 



3XiX3 {X2 - X3) {Xi - X2) 



Cfl(x) = exp , C^^{x) = exp 



3X3 {Xl - X2) {X3 - Xi) J ' \3Xi {X3 - X2) {X3 - Xi 



and all other C^Hx) as well as C{x) are equal to 1. 



6.8. Comparison with the elliptic Gamma cocycle ||F V2|| . In [[F V2|| , a nontrivial 



2-cocycle of SL(3,Z) with values in exp(27riQ(xj/xfc)[2;/x3]) was obtained by a similar 
construction involving the elliptic gamma function instead of U . In the language of 
discrete connections we use here, one defines an SL(3,Z) connection on the trivial line 
bundle over a dense subset of by setting 



X3 X3 X3 

Z X2- X3 X3 



^ I Z - X2 Xi - X2 Xi' ^ 

0r,3.2(a;, 2:) = T 

\Xi Xi Xi 

(l)r,i,jix,z) = 1, j^2. 

For our purpose, z may be considered here as a complex parameter on which SL(3, Z) 
acts trivially. The curvature of this connection was computed in ||F V2|| . It has the form 



^T,ij = 6xp(7riLfj(2;, x)), Cr(x) = 1, for some cubic polynomials L^j{x,z) in z with 
coefficients in Q(a:i, 0:2, X3). The relation to the curvature of (f) is 

Cij{x) = exp (-2m x Coefficient of in Lfj(a;, z)) . 

The curvature of a projectively flat connection defines an extension of the group and 
thus a characteristic class, see [ [b'V2|| . Conjecturally, C^j defines a nontrivial class in 
the group cohomology if^(SL(3, Z), exp(C(a;i, 0:2, X3))). This means (if the conjecture is 
true) that there is no gauge transformation given by exponentials of rational functions 
that can make the curvature trivial. 
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Appendix A. Theta functions and elliptic gamma functions 

We summarize some formulae about theta functions, gamma functions and phase 
functions, see |[1^'V2|| for more details. 

A.l. The theta function. Jacobi's first theta function is defined by the series 
e{z,T) = _^e*-a+i/2)^+2-*(^'+i/2)(^+V2), z,TeC, Imr > 0. 

It is an entire holomorphic odd function such that 

(27) 9{z + n + mr, r) = (^_iy'n+n^-mm^T-2.imzg^^^ ^ m, n G Z, 

and obeys the heat equation 

Its transformation properties with respect to SL(2, Z) are described in terms of genera- 
tors by the identities: 

0(^-z^r) = -0{z,t), 9{z,t + 1) = e^e{z,T), 6 (^^,-^^ = tV^e-^e{z,T). 

We choose the square root in the right half plane. 

A. 2. Infinite products. Let x,q EC with |g| < 1. The function 

oo 

{x;q) = Y[{l~xq^) 

3=0 

is a solution of the functional equation 

[qx-^q) = —^{x- q). 
1 — X 

Let X = e^""'' and q = e^^'^ . Then 

(28) e{z,T) = ie^'^^l^~'\x-q){q/x-q){q-q). 
We will also need the following variant of the theta function 6: 

eo{z,T) = {x; q){q/x; q) = -i— — ^(z,r). 

This function obeys 

eo{z + i,T) = eo{z,T), 

(29) eo{z + T,T) = -e~^^"eoiz,T), 

9o{t-z,t) = 6o{z,t). 

Its modular properties are Oq{z, r + 1) = 6q{z, r), and if z' = z/t, t' = — 1/t, 
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A. 3. Elliptic gamma functions. Here we consider two parameters r and p in the 
upper half plane, and set q = e^'^*'^, r = e'^'^^'^, and consider the function of a; = e^'^*^, 



(x; q, r) = Y\ (1 ~ xq^r^) = (x; r, q). 

j,k=0 



It is a solution of the functional equations 



{x]q,r) {x;q,r) 
{qx]q,r) = — — , {rx;q,r) 



(x;r) ' ' ' {x]q) ' 

The elliptic gamma function ^ , ||FV2|| is 

r 2;,r,p = — — = T{z,p,T). 

[x; q,r) 

It obeys the identities 

(30) T{z + l,T,p) = T{z,T,p), T{z + p,T,p) = eo{z,T)T{z,T,p), 

and is normalized by r((r + p + l)/2, r, p) = 1. The zeros of {x; q, r) are at x = q~^r~^, 
j,k = 0,1,2,.... They are all simple. Thus F has only simple zeros and simple poles. 
The zeros are at z = (j + l)r + (fc + l)p + 1, and the poles are at z = —jr — kp + l. Here 
j, k run over nonnegative integers and / over all integers. 

In fact, T{z, T,p) is, up to normalization, the unique 1-periodic meromorphic solution 
of u{z + p) = 6o{z, t)u{z) holomorphic in the upper half plane. 

A. 4. Modular properties. We consider the transformation properties of the elliptic 
gamma function under modular transformations of p and r. We have the identities 

T{z,T,p) = T{z,p,t), 
T{z,T + l,p) = T{z,T,p), 



T{z,T + p,p) 



T{z,T,p) 



T{z + T,T,p + t) 

r{z/p,T/p,-l/p) = e'-Q^^-'^'P^r{{z-p)/r,-l/r,-p/r)r{z,T,p), 

^, , z'^ r + p-1 2 t"^ + p^ + 'irp - 2,T - 'ip + I 

Q{z;t,p) = z + z 

Srp Irp brp 

+ ^{r + p-l)iT-'+p-'-l). 

A. 5. Extending the range of parameters. Since many operations we perform do 
not preserve the upper half plane, it is important to extend the range of values r and p 
can take. We set 

(x; q'^) = ^ (x; q-\ r) = - — - — -, (x; q, r'^^ - 



{qx;q) {qx;q,r) {rx;q,r) 
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These formulae define an extension of the functions {x; q), {x; q, r) to meromorphic func- 
tions on {(x, Q', r)||Q'| 7^ 1 7^ \r\}. It is clear that the functional relations 

{qx,q) = -^—{x,q), {qx]q,r) = -^—{x;q,r), 

still hold in this larger domain. Correspondingly, we extend the definition of and the 
elliptic gamma function by using the same formulae in terms of the infinite products. 
We obtain: 

111 

9o{z,-t)^— ■ r, r{z,-T,p)^— ■ r, r(z, T, -p) = 



0o{z + T,Ty ' ' '^^ r{z + T,T,py ^ ' ' r{z+p,T,py 

An easy check gives the following result: 

Proposition A.l. All identities for T and Oq of the preceding subsections continue to 
hold for all z, r, p such that r, p R. 

However, the statements about the position of zeros and poles are no longer valid. 

A. 6. The phase function. We keep the notation of the previous subsection and intro- 
duce a new variable a, and set a = e^'^*". The phase function is 

^ T{z + a,T,p) ^ {qr/xa;q,r){x/a;q,r) 

" ' ' r{z — a,T,p) {xa; q,r){qra/x; q,r) 

We have 

(31) n^{z + p,T,p) = + ^\ n^{z, T,p) ^ e'^ l^' + ^\ n^{z, t,p). 

uq[z — a,T) u[z — a,T) 

The properties of this function follow from those of the gamma function: 

Proposition A. 2. The function fla{z,T,p) obeys the identities 

(32) n,{z+p,T,p)^e'-^^^Jf±^Mz,T,p), 

(33) n,(z + T,T,p) = e'^^^^Jf±^Mz,T,p), 

U[z — a,p) 

(34) Q„(2;-M,T,p) = Qa(2;,T,p), 

(35) f^„(^,T,p) = 0„(^,p,r), 

(36) Qa{z,T,p) = Qa{z,T,T+p)Qa{z+P,T + p,p), 

(37) na{z,T,p) = riaiz + T,T,T + p)na{z,T + p,p), 

(38) Qa{z,T + l,p) = Qa{z,T,p+l) = Qa{z,T,p), 
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\T T T / 



Sa{z;T,p) = —{Qz^-6{T+p-l)z+2a^+T^+p^+3Tp-3T-3p+l), 

oTp 

(40) na{z,T,p)n_a{z,T,p) = l, 

(41) ,,,)|li±i^||i±^. 

t'^z — a, T)U{z — a,p) 

Appendix B. Modular properties of _R-matrices and qKZB operators 

Here we summarize some formulae giving the transformation properties of the R- 
matrix. Fix Ai, A2 G C and let R{zi — Z2, X, r, rj) be the i?-matrix of E-r^nish) associated 
to the evaluation Verma modules Vai(-2i), Vi^^{z2)- We have 

R{zi - Z2, A, r + 1, 77) = R{zi - z^, A, r, 77), 

and 

fl(£.^._i,!l) = e«-«'^^(A-2„;,(=))A,(A) 

\T T T T J 

X R{zi - Z2, A, r, r/)Ai(A)-M2(A - 2r//i«)-i. 
Here Ai{X) = A{zi, h^'\Ai, A, rj) with 



^7r 

^4(2;, /i, A, A, 1]) = exp — 

T 



z{hX - r^/i^) + ^(A - h){X + 7]K- r]h){X - r]A - r]h) 



These formulae can be deduced from the functional realization of representations ||FTV1 
the i?-matrix may be defined as the unique linear map such that 

R{zi - Z2, A, r, ri)Lj{zi, Z2, A, r, r]) = uj'^{zi, Z2, A, r, r]). 

The weight function uj = uiijef^ ® ef^ and the mirror weight function 0;^ take values in 
the tensor product Va^ ® Vx^ of evaluation Verma modules. They are given explicitly in 
terms of ratios of theta functions, see ||FTV1| , |FV1|] . 

From the above formulae, we deduce the transformation properties of the qKZB op- 
erators. One finds 

f^'^'-^'-'^) = e'^^U^^^^(--^^^^'^^^flAiz,-p9^^^ 

n j-1 

X Ki{z, A, T,p, rj) H A{z,-pei„ h^^\Aj, X - 2r]J2 v)'', 
j=i 1=1 
where we set 6ij = 1 if i > j and 6ij = ii i < j. In this somewhat simplified notation, 
the operators A depend on A and are to be viewed as multiplication operator by these 
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functions of A; Ki taken at A/r really means that we conjugate Ki by the operator of 
dilation of A by a factor 1/r. 

Appendix C. Estimates for theta functions 

Lemma C.l. (i) For all t in the upper half plane, there exists a constant Ci(r) > 
such that for all A G C, 



|^(A,r)| <Ci(r)exp TT 



Imr 



(ii) For all e > and r in the upper half plane, there exists a constant C2{t) > such 
that for all A such that min{ | A — r — rs | , r, s G Z} > e, 

\9{\,T)\>C,{T)eexp (vr^^^ 

Proof: Let S{t) be the compact set S{t) = {A G C | |Re(A)| < 1/2, |lmA| < lmr/2}. 
Then ^(A, r) vanishes in S{t) only for A = 0. 

We use the functional relation (^) of the theta function. If A = Aq + m' + mr, with 
m', m G Z and Aq G S{t) then 

Im A — Im An /t \ , t a \ 

= ° (i"^^°+i"^^)|g(Ao,r)| 
= e"^e-"^^|e(Ao,r)|. 

Then the claim follows by setting Ci(r) = maxje"'^* im?' |6'(Ao, r)|, Aq G 5'(r)} and 
C2{t) = min{e^'^* Im?' |^^(Ao, t)/Ao|, Aq G ^(r)}, so that when |Ao| > e, we have 

ei\o,r)>9i\o,r)j^>C2iT)e. 

1^0 1 

□ 
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